One-dimensional Anderson Localization: Devil's Staircase of Statistical Anomalies. 
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The statistics of wavefunctions in the one-dimensional (Id) Anderson model of localization is 
considered. It is shown that at any energy that corresponds to a rational filling factor / = - there 
is a statistical anomaly which is seen in expansion of the generating function (GF) to the order 
q — 2 in the disorder parameter. We study in detail the principle anomaly at / = | that appears in 
the leading order. The transfer-matrix equation of the Fokker-Planck type with a two-dimensional 
internal space is derived for GF. It is shown that the zero-mode variant of this equation is integrable 
and a solution for the generating function is found in the thermodynamic limit. 

PACS numbers: 72.15.Rn, 72.70.+m, 72.20.Ht, 73.23.-b 



-Introduction. Anderson localization (AL) enjoys 
an unusual fate of being a subject of advanced re- 
search during a half of century. The seminal paper by 
P.W.Anderson 1 opened up a direction of research on the 
interplay of quantum mechanics and disorder which is of 
fundamental interest up to now 2 . The one-dimensional 
tight-binding model with diagonal disorder -the Ander- 
son model (AM)- which is the simplest and the most 
studied model of this type, became a paradigm of AL: 
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In this model the hopping integral is deterministic £j = t 
and the on-site energy 54 is a random Gaussian vari- 
able uncorrelated at different sites and characterized by 
the variance ((5ei) 2 ) — w. The dimensionless parameter 
a 2 = w/t 2 describes the strength of disorder. 

The best studied is the continuous limit of this model 
in which the lattice constant a — > at ta 2 remaining 
finite^^. There was also a great deal of activity^ 
aimed at a rigorous mathematical description of Id AL. 
However, despite considerable efforts invested, some sub- 
tle issues concerning Id AM still remain unsolved. One of 
them is the effects of commensurability between the de- 
Broglie wavelength A e (which depends on the energy E) 
and the lattice constant a. The parameter that controls 
the commensurability effects is the filling factor / = j^- 
(fraction of states below the energy E). 

It was known for quite a while^ that the Lyapunov 
exponent (which is essentially the inverse localization 
length) takes anomalous values at the filling factors equal 
to i and i (compared to those at filling factors / be- 
yond the window of the size a 2 < 1 around f = h and 
/ = i, see Fig.l). Recentlj*^^ it was found that the 
statistics of conductance in Id AM is anomalous at the 
center of the band that corresponds to the filling factor 
/ = |. We want to stress that all these anomalies were 
observed for the AM Eq. (TTJ) in which the on-site energy 
is random. This Hamiltonian does not possess the chiral 
symmetry 2 -^ which is behind the statistical anomalies at 
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FIG. 1: Schematic representation of statistical anomalies in 
the localization radius. The dashed line represents the "bare" 
localization radius £0 = a^- sin 2 (71-/); the circles give the 
localization length at a filling factor being a simple fraction 
and solid lines give a sketch of behavior in the perturbed 
window. 



the center of the band E — in the Lifshitz model de- 
scribed by Eq.((T]) with the deterministic = and a 
random hopping integral tj. Thus the statistical anoma- 
lies in the Id AM raise a question about hidden symme- 
tries which do not merely reduce to the two-sublattice 
divisio n 2 ' 11 ^ 2 . Finally, the sign of the anomaly is dif- 



ferent for the center of the band f = h and the filling 
factor / = |. All these observations point out to a new 
phenomenon of the devil's staircase type which is essen- 
tially due to disorder and its interplay with the Bragg 
scattering off the underlying lattice. 



In this Letter we study the generating function (GF) 
$(m, 4>; x) that allows to compute all local statistical 
properties of Id AM. The simplest of them is the statis- 
tics of the wavefunction amplitude |^(x)| 2 characterized 
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by the moments I m = (\t/j(x)\ 2m ) 
2 
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where L is the total length of the system and to = 
a— sin 2 (tt/) is the "bare" localization length. 

We derive the corresponding transfer-matrix equation 
(TME) for the GF 



d x $(u, (j>; x) = [Lf(u, 4>) - u] $(it, 0; x), 



(3) 



which zero-mode variant </>; x) = </>) is indepen- 
dent of the space coordinate x) appears to be a partial 
differential equation (PDE) depending on two variables. 
One of them (denoted by u)is associated with the ampli- 
tude of the wavefunction tp ~ ^/ucos^ while the other 
(denoted by </>)has a physical meaning of its phase. The 
well known- TME in the continuous limit / <1 can be 
obtained by the averaging of this PDE over the phase 
variable <j) thus reducing it to an ODE in the single vari- 
able u. 

We show that there are statistical anomalies at any ra- 
tional filling factor / = 2. Namely, the operator Lf(u, 4>) 

in Eq.Q expanded in the disorder parameter a 2 



L f {u,4>) 



a 2 L { p{u,(j>) 



is such that 
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(4) 
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contains a regular part L^™' rog ' ) with a smooth depen- 
dence on / and an anomalous part that appears only at 

/ = 2Tn' 2Tn> "*§+n - In th(3 leadin S 0rder (" = °) in 

a 2 the anomalous term appears only at / = i. In the 
next order one can observe anomalies at / — ^ and g, 
etc. Though anomalous terms corresponding to the de- 
nominator q > 2 are small at weak disorder, they have 
an abrupt dependence on /. This allows to speak about 
the "devil's staircase of anomalies". 

We study in detail the principal anomaly at / = h. 
Remarkably, the corresponding zero-mode TME appears 
to be integrable. We find a unique solution to this equa- 
tion which describes any local statistics of wavefunctions 
in the center of the band. 

- Derivation of the TM equation. The starting point 
of our analysis is the TM equation for the generating 
function &j(u, 4>) on the lattice site j: 

<f> j+1 {u,cf>) = + ^ [£(u,(f))-ci((f>)u] \ $j(u,0-7r/), 

(6) 

where C(u,(f>) = C2{4>)u 2 d 2 + C3((j>) (ud u — 1) + 
Ci{4>) ududtj, + cs(</>) dff, + C(i((j)) d^. The coefficients Ci(4>) 



are all combinations of cos(20) and sin(2</>) which at 
first glance do not show any nice structure: ci(</>) = 
|(1 + cos(2</>)), c 2 (<P) = 1 - cos 2 (20), c 3 (0) = -(1 - 
cos(20) - 2cos 2 (2</>)), a((/>) = sin(20)(l + cos(2^)), 
c 5 (0) = -| sin(20)(l + cos(20)), c 6 (</>) = ±(l + cos(20)) 2 . 

This equation has been derived in Ref? 1 ^ by expan- 
sion to the first order in a 2 of the exact integral TM 
equation obtained by the super-symmetry method". By 
construction 1 ^ the function &j(u, 4>) must be periodic in 
<j) with the period of 7r which corresponds to the phase 
factor cos <f> of the wave function sweeping all possible 
values in the interval [0,7r]. However, the shift in the 
argument 4> m the r.h.s. of Eq.([6]) is by a fraction f of 
7r. For a rational / = | one has to make q iterations in 
Eq.© in order to get a closed equation for the GF. In 
the leading order in a we obtain: 
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The reason for the anomaly is the following identity that 
shows a jump at q = 2: 



9-1 



9-1 



^^ e 2i(p~2irTrp/q _ q ( 



Aiip—Air-irp/q _ \ 0; q > 2 

" l2e 4 ^, q = 2 

r=0 r=0 

(8) 

One can see from this identity that for q > 2 the sum- 
mation in Eq.((7|) is the same as averaging over <\>: all 
the (^-dependent terms vanish in both cases. Assuming 
q <C ^o/ a i expanding the l.h.s. of Eq.© and introducing 
the dimcnsionless coordinate x — ja/£o we obtain: 
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(9) 



This equation admits the independent of </> stationary 
solution: 



$(w, (f) = e 



;*VT=4i(2\AZ). 



(10) 



This solution has been earlier obtained^ in the continu- 
ous limit /< 1. It also arises in the theory of a multi- 
channel disordered wir o 14 i 15 . For a system of the size 
L — > oo only zero mode solution corresponding to e = 
is relevant. Substituting Eq. (TlT)]) with e = into Eq.® 
we found the following distribution function of the eigen- 
function amplitude in a long strictly one-dimensional sys- 
tem (amazingly, this result was not known before): 



la exp( 
L 



(ii) 



This distribution is valid for I'f/'PA) ^ e and should 
be cut off at very small \4>\ 2 to ensure normalizabilityi^. 
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At q = 2 (and only at q = 2 in the leading order in a) 
the ^-dependence in Eqs.©,® survives and gives rise 
to the anomalous terroi^: 



AL (0) = cos(#) 
+ sin(4</>) 



-u 2 dl + 2ud u + \d 2 ~2 



(12) 



Again, like in Eq.®, there is apparently no nice struc- 
ture in Eq. (fT2"l) . Moreover, because of the anomalous 
term the entire operator I/\ (u,<f>) acquires an explicit 

2 

^-dependence and thus the zero-mode TME becomes a 
two- variable second-order PDE which no longer admits a 
^-independent solution. Yet it appears exactly solvable! 

-Separation of variables. The integrability of the zero 
mode TME for / = | is shown in three steps. The 
step one is to introduce new set of variables u and 
v = u cos(2</>) instead of (it, (ft) and a new function 
$(w, v) = u^ 1 $(u, \ arccos(w/w)). In these variables the 
stationary TME [lA°' res) + AZ,(°)]$ = -e$ takes a very 

2 

symmetric form: 
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(13) 



where the operators D\ and D 3 belong to the family of 
three operators from the representation of the sl 2 algebra: 



Di = \Ju 2 — v 2 d u , D 2 = ud v +vd u , D 3 = — \J u 2 — v 2 d v 

(14) 

obeying the commutation relations: 

[D 1 ,D 2 ] = -D 3 , [Dz,D x ]=D 2 , [D 2 ,Dz]=D x . (15) 

Now it is clear that there is a hidden order in a set of 
^-dependent terms in Eq. (|12p and the way they match 



the regular part Lj ' 10 ^ in r.h.s. of Eq.([9]). 

The next step is to transform Eq. (fT3]) to the 
Schroedinger-like equation — (<9 2 + d 2 ) ^ + U(u, v) = 
for the function ^(u, v) = (u 2 — v 2 )* $, where 



U 



3 u 2 + v 2 



1 u 



4 {u 2 — v 2 ) 2 2 u 2 — v 2 u 2 — v 2 
Finally we introduce the variables 



(16) 



u + v 



6 = — ^— = u cos 2 <f), r\— U ^ — u sin 2 cj). (17) 

It is easy to see that the kinetic energy and the first 
two terms in Eq. (ll6|) become the sum of two identical 
Hamiltonians + H n where H% is given by: 



-a. 
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16 C 2 + 4f ' 



(18) 



Thus at e = we explicitly separated the variables 
in Eq. (|T5|) reducing the two-dimensional PDE to two 



ODE's of the Schredinger type H^<f\(^) = \tp\(£) and 
H n <p-\{rj) = —\(p-\(r/). Each of this equations reduces 
to the well known Weber's differential equation which so- 
lution is given in terms of the hypergeometric functions 
(Whittaker functions)^. 

Albeit the above procedure does not work for e ^ 
because of the last term in Eq. (fll)]) . the integrability of 
the zero-mode TME is a remarkable fact that allows to 
describe anomalous statistics in an infinitely long system. 

-Uniqueness of the solution. The general solution to 
the " Schroedinger equation" with e = in Eq. (fT6|) is 
given by the integral over the parameter A: 



dXdX c(A, A) (p\(Q tp-\(rj), 



(19) 



where integration is generically over the complex plane 
of A and c(A, A) is an arbitrary functioniS,. How does this 
huge degeneracy comply with the intuitive expectation 
that the statistics of wavefunctions in an infinite disor- 
dered chain should be unique and independent of the 
boundary conditions? Below we show that the natural 
physical requirements on $(u, (f) help to determine GF 
up to a constant factor which can be further fixed using 
the wave function normalization (l^l 2 ) = \- 

First of all we note that F(X; £, rf) — <p\(£) <p-\{rj) is a 
holomorphic function of A, i.e. it depends only on A = 
pe %cr but not on A = per la ' . The idea is to represent the 
integral over the complex plane as an integral over p and 
a and then rotate the contour of integration p — > te~ la 
so that the dependence on a remains only in c(A, A) and 
in the integration measure but not in ^(A;^^). Then 
performing integration over a one obtains a new function 
C(t) = t ] do- e~ 2ia c{t 1 te- 2la ) which stands for c(A, A) 
in an expression similar to Eq.(|19jl but involving only a 
one-dimensional contour integral. This contour can be 
further rotated to make the expression more symmetric. 
Thus without loss of generality we write a solution to the 



(20) 



zero-mode TM equation Eq. lfTBj) for / = -: 



Ac, 



V 



(^)l/4 



dAC(A) 
V4A/ 



+ c.c 



Here W K ^(z) is the Whittaker function^; e = e™/ 4 , e = 
e _l7r / 4 , and C(A) is a real function yet to be determined. 

Before we proceed with determining this function it is 
important to establish its properties as A — > 0. To this 
end we note that the integral over u in Eq. ^ with m = 1 
must be divergent at small u. Indeed, if it is convergent 



then the factor 



makes the first mo- 



(m-2)! r(m-l) 

ment equal to zero which contradicts the normalizability 
of wave function J dx(\4>(x)\ 2 ) = 1. At the same time 
for m > 1 the integral should converge to ensure finite 
higher moments. This implies that $(m, </>) must tend to 
a constant as in O^S. Given the asymptotic behavior 
of Whittaker functions this is equivalent to: 



C(A) = A"? C(A), (7(0) = const. 



(21) 
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GF defined by Eq. (|20|) is periodic in cj) with the period 
^ as it should be for q = 2. This is guaranteed by the 
adding of the c.c term in Eq . ([20| . What is not automat- 
ically guaranteed is that rf) is smooth as a function 
of (j> at <f) = 0. We will see that it is the requirement of 
smoothness at (j) — which fixes (up to a constant factor) 
the unknown function C(A). 

Indeed, the discontinuity of derivatives at <f> — may 
arise from the branching of the expression in Eq. ([20]) at 
a small ?y. From the representation of the Whittaker 
function in terms of the hypergeometric functions one 
concludes that the general solution Eq. (|2TJ)) is a sum of a 
part which is regular in the vicinity of r) = and a part 
which has a square-root singularity yjr\ ss -y/?I|</>|. The 
condition that this latter part cancels out in the solution 
Eq.® is the following (i is real): 



C(et) _ m /3 .3 i V 

TT —T- e 8t iF\ l it,-, — 

T(\-it) V 4 2'4f 



0.(22) 



The crucial fact for the possibility to fulfil this condition 
is the identity for the hypergeometric functions^: 



(23) 



A 1 2' 

Now one can immediately guess the solution for C(A): 

(24) 

It is easy to see that the general solution to Eq. (f!?2|) is 

00 

C(X)=C (X)S(X)=C (X) ^T^A 4 *, (25) 



(UX) = l ( \+e\\ rQ + eA 



fc=0 



where the function 5(A) must be regular in the entire 
complex plane of A. Now we apply the condition of con- 



vergence of the integral over A in Eq. (|20| at large A to find 
the allowed asymptotic behavior of 5(A) at A — > 00. Sub- 
stituting Eq. p5|) into Eq. (|20|) and using the asymptotics 
of the Whittaker and T-functions we find that the inte- 
grand behaves as A~ 3 5(A) at A — > 00. This means that 
1 5(A) I should increase not faster than A 2 . There is only 
one such entire function with the structure of Eq. ([25| : 
this is a constant 5(A) = oq = const. 

-Conclusion and discussion. Eqs. (|2T)| . (f2"l"j) . (f2"4"|) is the 
main result of the paper. They give an exact and unique 
(up to a constant factor) solution for the generating func- 
tion at / = \ anomaly which can be used to compute all 
local statistics of the one-dimensional Anderson model at 
L ->• 00. The integrability of TME Eq.© suggests that 
there is a hidden symmetry of the problem at / = | . We 
make a conjecture that this symmetry is naturally formu- 
lated in the three dimensional space rather than in the 
two-dimensional space (£, 77) and that it has to do with 
the symmetry of the 3d harmonic oscillator. This con- 
jecture is based on an analogy between our main result 
Eci. (f20|) and the expression for the Green's function of the 
3d harmonic oscillator problem 21 . This analogy concerns 
the parameter (A in our problem and k in Ref. 21 ) entering 
both in the argument of the Whittaker functions and in 
its first index in a mutually reciprocal way, as well as the 
second index of the Whittaker functions being i in both 
cases. Establishing this symmetry would also be useful 
for studying the anomalies at / = ~ with q > 2. We have 

obtained^! the operator iSP (u, (f>) in Eq.Q and shown 
that the mechanism similar to Eq.© is responsible for 
the anomaly at / = |. The results of this study will be 
published elsewhere. 
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